Spin-Hall effect on edge magnetization and electric conductance of a 2D 
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The intrinsic spin-Hall effect on spin accumulation and electric conductance in a diffusive regime 
of a 2D electron gas has been studied for a 2D strip of a finite width. It is shown that the spin 
polarization near the flanks of the strip, as well as the electric current in the longitudinal direction 
exhibit damped oscillations as a function of the width and strength of the Dresselhaus spin-orbit 
interaction. Cubic terms of this interaction are crucial for spin accumulation near the edges. As 
expected, no effect on the spin accumulation and electric conductance have been found in case of 
Rashba spin-orbit interaction. 
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Spintronics is a fast developing area to use electron 
spin degrees of freedom in electronic devices . One of 
its most challenging goals is to find a method for manip- 
ulating electron spins by electric fields. The spin-orbit 
interaction (SOI), which couples the electron momentum 
and spin, can be a mediator between the charge and spin 
degrees of freedom. Such a coupling gives rise to the so 
called spin-Hall effect (SHE) which attracted much in- 
terest recently. Due to SOI the spin flow can be induced 
perpendicular to the DC electric field, as has been pre- 
dicted for systems containing spin-orbit impurity scat- 
terers 0. Later, similar phenomenon was predicted for 
noncentrosymmetric semiconductors with spin split elec- 
tron and hole energy bands |3l ■ It was called the intrinsic 
spin-Hall effect, in contrast to the extrinsic impurity in- 
duced effect, because in the former case it originates from 
the electronic band structure of a semiconductor sample. 
Since the spin current carries the spin polarization, one 
would expect a buildup of the spin density near the sam- 
ple boundaries. In fact, this accumulated polarization is 
a first signature of SHE which has been detected experi- 
mentally, confirming thus the extrinsic SHE j3| in semi- 
conductor films and intrinsic SHE in a 2D hole gas j^. 
On the other hand, there were still no experimental evi- 
dence of intrinsic SHE in 2D electron gases. The possibil- 
ity of such an effect in macroscopic samples with a finite 
elastic mean free path of electrons caused recentl y rn uch 
debates. It has been shown analytically 
and numerically 12] that in such systems SHE vanishes 
at arbitrary weak disorder in DC limit, for isotropic, as 
well as anisotropic ^3 impurity scattering, when SOI is 
represented by the so called Rashba interaction • As 
one can expect in this case, there is no spin accumulation 
at the sample boundaries, except for the pockets near the 
electric contacts "T]. At the same time, the Dresselhaus 
SOI 01, which dominates in symmetric quantum wells, 
gives a finite spin-Hall conductivity • The latter can 
be of the order of its universal value e /SttH. The same has 
been shown for the cubic Rashba interaction in hole sys- 



tems [l^ . In this connection an important question 
is what sort of the spin accumulation could Dresselhaus 
SOI induce near sample boundaries. Another problem 
which, as far as we know, was not discussed in literature, 
is how the electric current along the applied electric field 
will change under SHE. In the present work we will use 
the diffusion approximation for the electron transport to 
derive the drift diffusion equations with corresponding 
boundary conditions for the spin and charge densities 
coupled to each other via SOI of general form. Then the 
spin density near the flanks of an infinite 2D strip and 
the correction to its longitudinal electric resistance will 
be calculated for Dresselhaus and Rashba SOI. 

Let us consider 2DEG confined in an infinite 2D strip. 
The boundaries of the strip are at y = ±d/2. The elec- 
tric field E drives the DC current in the 2:-direction and 
induces the spin-Hall current in the y-direction. This 
current leads to spin polarization buildup near bound- 
aries. Since d ^ kp^, where kp is the Fermi wavector, 
this problem can be treated within the semiclassical ap- 
proximation. Moreover, we will assume that d is much 
larger than the electron elastic mean free path I, so that 
the drift-diffusion equation can be applied for description 
of the spin and charge transport. Our goal is to derive 
this equation for SOI of general form 



Hso hk ■ (T , 



(1) 



where (T={a^ , ,a^) is the Pauli matrix vector, and the 
effective magnetic field hk = h k is a function of the 
two-dimensional wave- vector k. 

We start from determining linear responses to the mag- 
netic B{r,t) and electric V{r,t) potentials. The mag- 
netic potentials are introduced in order to derive the dif- 
fusion equation and play an auxiliary role. The corre- 
sponding one-particle interaction with the spin density 
is defined as Hgp — B{r,t) ■ a. These potentials induce 
the spin and charge densities, S{r,t) and n{r,t), respec- 
tively. Due to SOI the charge and spin degrees of freedom 
are coupled, so that the electric potential can induce the 
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spin density and vice versa. Therefore, it is conve- 
nient to introduce the four- vector of densities Di{r,t), 
such that Do{r) = n{r,t) and D.^^y^^{r,t) = S^^y^;,{r,t). 
The corresponding four-vector of potentials will be de- 
noted as ^i{r,t). Accordingly, the linear response equa- 
tions can be written in the form 



(2) 



The response functions Ily (r , r', t — t') can be expressed 
in a standard way |17i] through the retarded and ad- 
vanced Green functions G'^{r,r',t) and G'^{r,r',t). In 
the Fourier representation we get 

xE,G'^(r,r',c^' + c^)E,]), (3) 

where Sq = 1, = (Ji at z = x,y,z and npiuj) is the 
Fermi distribution function. The time Fourier compo- 
nents of densities Df{r, t) at w ^ Ep are defined as 

J ^ J 

X (Tr[G'-(r,r',c^')S^G''(r',r,c^')S, 

- Gnr,r^c^')S^Gnr^fc,c^')S,]). (4) 

The trace in Eqs. (|3I4|I runs through the spin variables, 
and the angular brackets denote the average over the ran- 
dom distribution of impurities. Within the semiclassical 
approximation the average of the product of Green func- 
tions can be calculated perturbatively. Ignoring the weak 
localization effects, the perturbation expansion consists 
of the so called ladder series [itIIi^ . At small uj and large 
r — r'\ they describe the particle and spin diffusion pro- 
cesses. The building blocks for the perturbation expan- 
sion are the average Green functions and Q"" , together 
with the pair correlator of the impurity scattering poten- 
tial Uscir). A simple model of the short-range isotropic 
potential gives {Usc{r)Usc{i'')) = T6{i — r')/TTNo, where 
Nq is the electron density of states at the Fermi energy 
and r = l/2r. Within the semiclassical approach the 
explicit behavior of the electron wave functions near the 
boundaries of the strip is not important. Therefore, the 
bulk expressions can be used for the average Green func- 
tions. Hence, in the plane wave representation 

g^{k, Lo) = [g'^ik, Lo)]^ = {uj-Ek-hk-c7 + iT)-\ (5) 

where Ek—k^/{2m*) — Ep. Since the integral in Q 
rapidly converges at \r — r'\ < kp^, D^{r,u)) are given 
by the local values of potentials. From Q-© one easily 
obtains the local equilibrium densities 



In their turn, the nonequlibrium spin and charge densities 
are represented by the first term in Eq. Within the 
diffusion approximation this term is given by the gradient 
expansion of Q [l^ . Such an expansion is valid as far as 
spatial variations of Di{r,u;) are relatively small within 
the length of the order of the mean free path I. The 
length scale for spin density variations near the bound- 
aries of the strip is given by vp/hkp. Hence, the dif- 
fusion approximation can be employed only in the dirty 
limit hkp <C l/r. The diffusion equation is obtained 
after the ladder summation in the first term of Eq. ^ 
and multiplying this equation by the operator inverse to 
Hij(r,r',w), as it has been previously done in [l^ l2fij|. 
After some algebraic manipulations one gets 



(7) 



where the diffusion operator 7)^^ can be written as 

^ S'WV^ - + i?'^"V™ + AP' . (8) 

The first term represents the usual diffusion of the spin 
and charge densities, while the second one describes the 
D'akonov-Perel' 0| spin relaxation 



(9) 



where i,j ^ 0, the overline denotes the average over the 
Fermi surface and = hk/hk- The third term gives 
rise to precession of the inhomogeneous spin polarization 
in the effective field of SOI ^ 



4t J2 hlv 



(10) 



The nondiagonal elements of the form I)'" appear due to 
spin-orbit mixing of spin and charge degrees of freedom. 
They are collected in M'^ . For Rashba SOI M'° have 
been calculated in 0, ■ In general case 



hi dnl. 
r2 dk 



(11) 



D^{r,Lo) = -2No^^{r,uJ). 



(6) 



When a time independent homogeneous electric field 
is applied to the system one has $o = sEx and Dq — 
~2NQeEx. At the same time, = and, hence, ^ 
at i = x,y,z. Due to charge neutrality the induced 
charge density cDq = 0. It should be noted that in the 
system under consideration the charge neutrality can not 
be fulfilled precisely. The spin polarization accumulated 
at the strip boundaries gives rise to charge accumulation 
via the terms in 0-®- The screening effect will, 
however, strongly reduce this additional charge, because 
the screening length of 2DEG is much less than the typi- 
cal length scale of spin density variations. We will ignore 
such a small correction and set Dq = in (TJ. In this 
way one arrives to the closed diffusion equation for three 
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components of the spin density. This equation coincides 
with the usual equation describing diffusive propagation 
of the spin density , for exception of the additional 
term -M'^D'^ = 2NoeEh^^l^/V'^ due to the external 
electric field. Its origin becomes more clear in an infi- 
nite system where the spin density is constant in space 
and only F*^ and M*-' are retained in JZ!)-©. Hence, the 
solution of © at = is 5'i = , with 



5' ^ A°/2 



Y2 l^y^ ) ^kQj^^ 



(12) 



where (r~^)'^ is the matrix inverse to Such a phe- 
nomenon of spin orientation by the electric field was pre- 
dicted in Ref. jl^] and recently observed in In the 
special case of Rashba SOI hk = a(k) (k x z) it is easily 
to get from (O the resuh of Rcf. 16] S*^ = -NoeEar. 

In addition to the diffusion equation one needs the 
boundary conditions. These conditions are that the three 
components of the spin flux 1^,1^,1^ flowing in the y 
direction turn to at y = id/2. The linear response 
theory, similar to l(2Jl gives 

4(r,t) = Jd\'dt'Y,Ei^{r,r',t~t')^,{r',t'), (13) 

j 

where the response function S is given by 
duj' dnF{uj') 



S^^.(r,r',c.) = tujj —^l^{Tr[G'^{r',r,u;') x 

xJlG-{r,r\u;' + u;)j:,]), (14) 

with the one-particle spin-current operator defined by 
J^ = ((tV;-|-v/(T*)/4 and the particle velocity 



(15) 



Taking into account (|7I5|I . we obtain from p3ll4|l 
.dS, 1 



Ifir) 



-D- 



dy 



R}^y{Sj-S^) + 5,JsH- (16) 



The first two terms represent the diffusion spin current 
and the current associated with the spin precession. The 
third term is the uniform spin-Hall current polarized 
along the z axis. It is given by 



IsH = -\R^-yS) + eE^vl{'^.h,).. (17) 



From (|10I12II it is easily to see that for Rashba SOI both 
terms in (I17II cancel each other making Ish = 0, in accor- 
dance with li HHP, To", Til [13. Therefore, in case of the 
strip the solution of the diffusion equation satisfying the 
boundary condition is Sj = SjySy. Hence, the spin den- 
sity is uniform and does not accumulate near boundaries. 
It should be noted that such accumulation can, however, 




FIG. 1; Spin densities AS'i(±d/2) = ASi for i ^ x,z on 
the boundaries of the strip, as functions of its width d, for 
K/fc=0.9,1.0, and 1.3, respectively. The insert shows the de- 
pendence of ASz{y) on the transverse coordinate y. Lengths 
are measured in units of Iso ~ Vph/2{vFyhhy) 



take place in the ballistic regime of electron scattering 
[2^. At the same time, as shown in Ref. even in 
the diffusive regime Ish / for the Dresselhaus SOI. 
This inevitably leads to the spin accumulation. Taking 
Dresselhaus SOI in the form 



/3kx {ky 



K 



Pky{k'^ 



(18) 



one can see that the bulk spin polarization (|12() has a 
nonzero 5*^ component, i?^^^ 7^ 0, while R^^y = 0. 
Hence, the solution of the diffusion equation with 
the boundary condition iy{±d/2) = /|'(±d/2) is 
5^, S;, ^0 Sy ^0. Let us define AS,{y) = S,{y) - Sf. 
The dependence of ASi{ztd/2) from the strip width, 
as well as an example of A^z coordinate dependence, 
are shown in Fig. 1. The damped oscillation in the d- 
dependence of the spin accumulation on the flanks of the 
strip can be seen for the Sz polarization. Similar oscilla- 
tions take place also in the coordinate dependence. The 
length scale of these oscillations is determined by the spin 
precession in the effective spin-orbit field. 

The arbitrary units have been used in Fig. 1. For a nu- 



merical evaluation let us take E = lO^V/m, 



r/h = 



0.1, and K/kp—0.8 for a GaAs quantum well of the width 
ui^lOOAdoped with 1.5 • lO^^m"^ electrons. We thus ob- 
tain \ASz{±d/2)\ ~ 5 X lO^^m^^. The corresponding 
volume density ASz/w ~ 5 x lO^^m"'^, which is within 
the sensitivity range of the Faraday rotation method 0] . 
It should be noted that in the considered here "dirty" 

limit \J'h^T/h <Si 1 the spin-Hall current is suppressed 
by the impurity scattering. As shown in |llLll2| for Dres- 
selhaus and cubic Rashba SOI, this current decreases as 
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/i|^r^/?i^ down from its highest universal value. At the 
same time, an analysis of the diffusion equation shows 
that the accumulated at the flanks of the strip spin den- 
sity decreases slower, as 



'h'^^r/h. This explains why 
for the considered above realistic numerical parameters, 
even in the dirty case, the noticeable spin polarization 
can be accumulated near the boundary. 

Usually, the spin-Hall effect is associated with the spin 
polarization flow, or the spin density accumulation on 
the sample edges, in response to the electric field. On 
the other hand, this effect can show up in the electric 
conductance as well. To see such an effect we take 0- 
projection of ifT^ . which by definition is the electric cur- 
rent. The current flows along the X axis. The correspond- 
ing response function is given by l|14|l with Jq = v^;. 
Using Eqs. ((T3|l . Q and expressing from Q one 

gets the electric current density 



r = (je + a 



dy 



where a is the Drude conductivity and 



2r2 
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dhk 



X h 



k)z 



(19) 



(20) 



The total current is obtained by integrating (|19|l over y. 
Therefore, the spin-Hall correction to the strip conduc- 
tance 

A 2 A 

AG = - {S,{d/2) - S,{-d/2)) = —S,{d/2) . (21) 

Hence, the dependence of AG on the strip width coin- 
cides with that of the spin density shown in Fig. la. 

In conclusion, we employed the diffusion approxima- 
tion to study the spin-Hall effect in an infinite 2D strip. 
In case of the Dresselhaus spin-orbit interaction this ef- 
fect leads to spin accumulation near the flanks of the 
strip, as well as to a correction to the longitudinal elec- 
tric conductance. Both, the spin accumulation and the 
conductance exhibit damped oscillations as a function of 
the strip width. 

This work was supported by the Taiwan National Sci- 
ence Council NSC93-2112-M-009-036 and RFBR Grant 
No 030217452. 



Samarth, eds., Semiconductor Spintronics and Quantum 
Computation (Springer- Verlag, Berlin, 2002); I. Zutic, 
J. Fabian, and S. Das Sarma, Rev. Mod. Phys. 76, 323 
(2004). 

[2] M. I. Dyakonov, V. I. Perel, Phys. Lett. A 35, 459 (1971); 

J. E. Hirsch, Phys. Rev. Lett 83, 1834 (1999). 
[3] S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, 

1348 (2003); J. Sinova et. al, Phys. Rev. Lett. 92, 126603 

(2004); D. Culcer et. al, Phys. Rev. Lett. 93, 046602 

(2004). 

[4] Y. K. Kato, et. al, Science 306, 1910 (2004) 
[5] J. Wunderlich et. al, Phys. Rev. Lett. 94, 047204 (2005) 
[6] J. I. Inoue, G. E. W. Bauer, and L.W. Molenkamp, 
Phys. Rev. B 70, 041303 (2004); E. I. Rashba, 
|cond-m at/0409476 (2004) (unpublished); O. Chalaev and 
D. Loss, cond-mat/0407342 (unpublished). 
[7] E. G. Mishchenko, A. V. Shytov, and B. I. Halperin, 

Phys. Rev. Lett. 93, 226602 (2004). 
[8] A. A. Burkov, A. S. Nunez, and A. H. MacDonald, Phys. 

Rev. B 70, 155308 (2004). 
[9] O. V. Dimitrova, cond-mat/0405339 (unpublished) 
[10] R. Raimondi, P. Schwab, Phys. Rev. B 71, 033311 (2005) 
[11] A. G. Malshukov and K. A. Chao, Phys. Rev. B 71, 

121308(R) (2005) 
[12] B. A. Bernevig, S.C. Zhang, co nd-mat/0411457j K. No- 
mura et al., cond-mat/0506189 
[13] Yu. A. Bychkov and E. I. Rashba, J. Phys. C 17, 6039 
(1984) 

[14] G. Dresselhaus, Phys. Rev. 100, 580 (1955) 

[15] The cubic Rashba interaction should not be confused 
with the conventional linear Rashba SOI with the wave- 
vector dependent coupling constant a{k). In the latter 
case SHE (xa{k)k/EF < 1 |3,|l2. 

[16] V. M. Edelstein, Solid State Commun., 73, 233 (1990); 
J. I. Inoue, G. E. W. Bauer, and L.W. Molenkamp, Phys. 
Rev. B 67, 033104 (2003) 

[17] A. A. Abrikosov, L.P. Gorkov, and I.E. Dzyaloshinskii, 
Methods of Quantum Field Theory in Statistical Physics, 
(Dover, New York, 1975) 

[18] B. L. Altshuler and A. G. Aronov, in Electron-Electron 
Interactions in Disordered Systems, edited by A. L. Efros 
and M. PoUak (North-Holland, Amsterdam, 1985). 

[19] A. G. Malshukov and K. A. Chao, Phys. Rev. B 61, 
R2413-R2416 (2000). 

[20] A. G. Mal'shukov, K. A. Chao, and M. Willander, Phys. 
Rev. Lett. 76, 3794 (1996); Physica Scripta T66, 138 
(1996). 

[21] M. I. D'yakonov and V. I. Perel', Sov. Phys. JETP 33, 
1053 (1971) [Zh. Eksp. Teor. Fiz. 60, 1954 (1971)]. 

[22] Y.K. Kato, et. al, condmat/0502627. 

[23] B. K. Nikolic, et. al, cond-mat/0412595 Q. Wang, L. 
Sheng, C. S. Ting, |Cond-mat /0505576| 



[1] G. A. Prinz, Science 282, 1660 (1998); S. A. Wolf et. al. 
Science 294, 1488 (2001); Awschalom, D. Loss, and N. 



